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Abstract 

This poster discusses a possible explanation for the relationship between the mass of the 
central supermassive black hole and the velocity dispersion in the bulge of the host galaxy. 
We suppose that the black hole and the dark matter halo are forming simultaneously as 
matter falls in, and a self-similar system then exists in which the mass and the velocities 
of the system evolve as power-law functions of time. This leads naturally to a relationship 
between the black hole mass and the velocities in the halo which, with a reasonable choice 
of cosmological parameters, is in good agreement with the observed relationship. We also 
confirm this relationship with more robust numerical results. 

1.1 Dynamical Growth of Black Holes 

Supermassive black holes (BHs) are now considered to be a common feature of 
galaxies which have a bulge. Furthermore, a number of observational properties of the host 
galaxy correlate with the BH mass. Among the strongest of these correlations is the rela- 
tionship between the BH mass and velocity dispersion within the galactic bulge (Ferrarese 
& Merritt 2000; Gebhardt et al. 2000): M BH oc a a , where a = 4.02 ± 0.32 (Tremaine et al. 
2002). Since the velocity dispersion is measured well outside of the BH "influence," this 
correlation indicates that an intimate relationship exists between the BH and the dynamical 
structure of the host galaxy. 

We assume that a galaxy forms by the extended collapse of a "halo" composed of col- 
lisionless matter and that simultaneously the central black hole is growing proportionally 
to the halo as matter continues to fall in. This is equivalent to the assumption of multidi- 
mensional self-similarity (Carter & Henriksen 1991; Henriksen 1997), and MacMillan & 
Henriksen (2002) show that, in this case, the mass M inside any surface in the halo is related 
to the velocity of the particles by 



where a is the averaged velocity. Since the system is self-similar, this relation will apply at 
both the "bulge" mass scale as well as the BH mass scale. The quantities 5 and a are scales 
in space and time, respectively, and are related to the power-law index of the initial density 
perturbation e in spherical infall models of halo growth (Henriksen & Widrow 1999): 
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This power-law index e depends on the spectral index n of the primordial power spec- 
trum, P(k) oc k", by n = 2e— 3. There is therefore a direct link between the initial power 
spectrum and the predicted relationship between the BH mass and the velocity dispersion; 
for agreement with observations, we must have n = -2, so that e = 1/2 and 5 /a = 2. 

We determine below that the addition of sufficient angular momentum breaks the self- 
similarity in the central regions. This creates a "mass reservoir" around the BH that does 
grow in proportion to the galaxy mass, and from which the BH grows more slowly by colli- 
sional interactions between clumps of matter. Provided most of this mass may be absorbed 
by the BH on a cosmic timescale, the proposed relation should still hold. 

1.2 Numerical Simulations 

As an extension to the analytical work of MacMillan & Henriksen (2002), we 
present here numerical simulations which confirm the predicted relationship between the 
BH mass and the velocity dispersion. 

The initial conditions for the simulations follow those of Henriksen & Widrow (1999): 

p( r ,td = P c(td[l+A(r,tdl (1-3) 

where 

^ tl) { A(l-B)(r/r c r r>r c . {lA) 

This initial density profile is that of a power-law function of radius, modified at small radii 
to account for suppression of small scale fluctuations. The constant B = 5e/(3e + 6) and we 
choose A = 0.5. The particles are started in the unperturbed Hubble flow. 

Our initial numerical work began with a simple shell code, which is used to follow spher- 
ically symmetric shells of matter in an Einstein - de Sitter universe; each shell evolves ac- 
cording to 

d 2 n _ GM(n) f 

where M(r,) is the total mass enclosed by the shell at radius r, . The angular momentum j 
is given by j 2 = 2JGM(ri)ri, evaluated at the initial shell positions. Thus the particles are 
allowed to have, in general, elliptical orbits, and the shells follow the evolution of their radial 
coordinates. 

To extend the spherically symmetric numerical results, we consider black hole growth in 
an n-body system of n ~ 2000 particles. We use a modified version of Piet Hut's code 
( http://www.ids. ias.edu^piet/act/comp/algorithms/starter/i This code was modified to in- 
clude a softening length to handle close encounters between particles. 

In both simulations, a central BH has been included that can "eat" particles which fall 
below the Schwarzschild radius, R s = 2GMbh/c 2 , and grow in mass accordingly. Also, 
we choose a power-law index for the initial perturbation of e = 1/2, so that 6/a = 2 and 
Mbh cx a 4 , as observed. 
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1.2.1 Results 

The shell code simulations are run, with an initial BH mass of 10" 4 , until t = 100, 
at which time all matter is essentially "eaten" by the BH in the 7 = case. Other values of 
the specific angular momentum are also explored. The BH mass, as a function of time, is 
shown in Figure[TTJi for J = 0, 0.01, 0.1, 0.2, and 0.5. 

MacMillan & Henriksen (2002) predict that the BH mass should grow as t A . This is indeed 
the asymptotic behavior for the case where 7 = 0, as shown in the figure. This logarithmic 
slope, however, decreases as the specific angular momentum increases; for 7 = 0.1, the mass 
is increasing as t 29 . For larger values of 7, the growth essentially "shuts off" at some time, 
indicating that shells with high initial angular momentum are never reaching the centre of 
the system. 

Three different cases are explored with the n-body simulations. The first is an initial 
spherically symmetric system. In the second and third cases, spherical symmetry is removed 
by giving each of the particles a random displacement, with amplitude Aq. Systems are 
considered in which the amplitudes are about 0.8% and 1.6% of the initial extent of the 
system. 

Results are shown in Figure II. lb for the black hole mass as a function of time. For 
the case with no initial random displacements (Ao = 0), the growth is approximately that 
expected, although at later times the black hole growth slows down due to a lack of infalling 
particles. For the more general cases where Ao > 0, black hole growth remains self-similar, 
but with a much lower logarithmic slope: Mbh oc f or a little greater. In these systems, the 
particles have enough angular momentum that they don't fall straight into the central black 
hole, and rather form a "core" at the centre which slowly feeds the black hole growth. 

To test the prediction that Mbh « a 4 , the velocities of some fraction of the particles are 
averaged (equating a = v). The predictions of MacMillan & Henriksen (2002) apply to a 
fixed scaled position vector X, implying that we should average all particles within some co- 
moving radius. Results are shown in Figure lT"2"h for a comoving radius of r = 0.25. Although 
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Fig. 1.2. The measured Mbh - o relation for (a) a co-moving radius, and (b) a fixed 
radius. 

the case where Aq = may have a regime where it follows the predicted Mbh _ & relation, 
the cases with Aq > clearly follow no relation at all; their average velocity is about 0.18 
for most black hole masses. 

Of interest, however, are the results if we take a fixed radius, rather than a comoving one, 
and average the velocities of the particles within that radius. The plot of Mbh versus a is 
shown in Figure fT~2b . All three cases clearly show a trend that goes like Mbh oc cr 4 - 

1.3 Discussion 

For systems in which there is no angular momentum or deviation from spherical 
symmetry initially, the numerical work, with both the spherical shell code and the more 
general n-body code, confirms the predicted black hole growth: Mbh oc ? 4 . However, for 
the simulations in which there was some deviation from a simple spherical infall of matter, 
results were contrary to what was expected. In particular, the n-body results show a black 
hole growth approximately proportional to the time. 

However, it is apparent from the simulations that the black hole growth in these cases 
is not dominated by the infalling matter, as it is in the initially spherical case. Rather, the 
particles form a "core" about the centre, and so the growth is fed in a different manner. 
Dimensional arguments suggest that the Schwarzschild radius should scale with time as 
R s oc ct , where c is the speed of light, a constant of the system. If this is the case, then we 
get that the black hole mass will grow as Mbh oc R s oc f, which is approximately what is 
observed. Note, however, that this breaks the symmetry of the system, since the mass on a 
larger scale is predicted to grow as f 35 " 2 ". 

The numerically derived Mbh -ct relationship, which here compared the black hole growth 
with the average velocity of particles within some radius r, only gave the expected relation, 
Mbh oc er 4 , when a fixed, rather than comoving radius was considered. Furthermore, the 
cases for which the growth rate was much shallower than expected still showed the same 
trend as the other case; that is, Mbh oc a 4 regardless of how the black hole grew with time. 

As stated above, however, this result is for a fixed radius. From the definition of X and T 
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in MacMillan & Henriksen (2002), we see that r = Xt S ' a . Fixing r requires X to then depend 
on time. Recalling that we're considering a density profile that goes as X~ e , and e is related 
to the quantity S/a by equation il.2\ . we can write the mass as a function of X: 



Thus, as a function of the (fixed) radius, we have M(r) oc r 3 2a / s . Similarly, we can write 
the self-similar velocity, assuming it is in its virialized state, as a function of X: 



so that the radial dependence of the velocity takes the form v(r) oc r l ~ a l s . 

Combining these two equations for mass and velocity and eliminating r gives us the fa- 
miliar relation given by equation 1 11. 11 1. Thus, taking the mass and velocity at a fixed radius 
r gives the same relation regardless of the time dependence of either. 

1.4 Conclusions 

We consider black hole growth in a system which forms as the black hole is grow- 
ing. If this is the case, then the system has multidimensional self-similarity, and the black 
hole should grow as a power-law function of time, Mbh oc t 3S ~ 2a . Furthermore, the veloci- 
ties of the particles of the system will also evolve in a similar manner, and we predict that 
log(M BH ) oc (3S/a-2)/(S/a- l)loga. 

Numerically, we confirm that the black hole does indeed grow as a power-law function of 
time; however, it only follows the predicted relationship if the initial system has little or no 
angular momentum. For systems with non-spherical initial conditions, the black hole grows 
as Mbh oc f or a little greater. 

Furthermore, we test theMBH-c relation with n-body simulations. Although results do 
not show a strong relation if we take an averaged velocity within a comoving radius, the 
predicted relation, that Mbh oc er 4 , is clearly present if we take a fixed radius and average the 
particles within that. This can be explained by considering the self-similar variable X to be 
a function of time such that the time dependence of the radial coordinate is fixed. 
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